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Hudson,  Colorado  so  that  data  could  be  collected  on  a large  two-dimension- 
al  array  in  a format  suitab’c  for  use  by  Raytheon  and  by  others  for 
adaptive  processing.  The  principal  modifications  included  the  construction 
of  a 32  element  cross  array  of  vertical  inonopoles  and  the  addition  of  quad- 
rature outputs  to  the  existing  32-channel  receiver  and  data  acquisition 
system.  Additional  Raytheon  tasks  included  the  development  of  software 
(for  controlling  the  data  collection  system,  for  processing  and  displaying 
the  field  data,  for  producing  simulated  field  tapes  for  specific  array 
geometries  and  model  environments),  data  collection  under  a variety  of 
conditions,  and  analysis  (numerical  simulations  to  guide  the  design  of  the 
antenna  array  and  data  reduction  to  evaluate  the1  performance  improvement 
obtained  through  adaptive  processing). 

The  principal  result  of  this  work  is  the  verification  that  the  HF  environ- 
ment will  in  fact  permit  most  of  the  performance  improvement  anticipated 
from  idealized  numerical  simulation,  of  adaptive  processing  applied  to 
large  two-dimensional  arrays.  Noise  floor  reductions  of  order  30  d 13  (com- 
pared to  conventional  processing)  were  realized  over  most  oi  the  radio  sky 
under  conditions  in  which  a single  strong  interferer  was  present  (10  dI3 
element  signal  to  noise  ratio).  Viev/ed  as  a technique  for  producing 
improvecl  two-dimensional  maps  of  the  radio  environment  (e.  g.  , radio 
location  or  noise  survey  applications)  rather  than  for  improving  signal  to 
noise  ratio,  the  maps  produced  by  adaptive-  beamforming  exhibit  superior 
angular  definition  and  apparent  freedom  from  sidelobe  responses.  That 
is,  the  adaptive  processor  is  effective  in  deconvolving  the  antenna  pattern 
from  the  measurements. 

The  HF  environment  is  found  to  change  rapidly  enough  so  that  if  a one  sec- 
ond block  of  data  is  processed  with  adaptive  weights  computed  for  the 
preceding  one  second  data  block,  approximately  10  dB  of  the  noise  floor 
improvement  attainable  by  using  optimum  weights  is  lost.  It  appears  that 
the  direct  matrix  inversion  method  for  implementing  adaptive  processing 
will  in  practice  require  a new  inversion  for  each  data  block  but  it  has  not 
yet  been  determined  how  long  these  data  blocks  can  be  made  without 
incurring  an  unacceptable  signal  to  noise  ratio  penalty. 

It  was  found  that  the  output  of  the  adaptive  processor  was  typically  20  dB 
less  than  that  of  the  conventional  processor  when  both  were  pointed  in  the 
source  direction,  an  unexpected  effect  (ideally  they  would  be  equal)  attri- 
buted to  the  departure  of  the  wavefront  of  the  incident  signal  from  an  ideal 
plane  wave.  This  departure  is  believed  due  to  ionospheric  multipath  but 
terrain  scattering  and/or  residual  calibration  errors  in  the  receiving  sys- 
tem cannot  be  ruled  out. 

The  result  of  the  array  geometry  study  is  to  suggest  low-redundancy  sparse 
arrays  as  promising  candidates  for  adaptive  antenna  applications  in  which 
both  angular  resolution  and  signal  to  noise  ratio  arc  considerations. 

It  is  recommended  that  the  processing  of  the  data  base  already  on  tape  be 
completed  so  as  to  explore  the  impact  of  variable  data  block  length  on  sig- 
nal to  noise  ratio,  to  explore  the  consequences  of  simultaneous  adaptivity 
in  azimuth,  elevation,  and  frequency,  and  to  explore  the  performance  of 
the  adaptive  processor  in  complex  interference  environments. 
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ABSTRACT 


The  DARPA  HF  Adaptive  Array  Program  has  been  a one  year,  three  com- 
ponent effort  (University  of  Colorado,  Stanford  Research  Institute,  and  Raytheon) 
to  determine  the  applicability  of  known  adaptive  antenna  processing  techniques 
to  HF  systems.  This  report  is  the  final  report  on  the  Raytheon  portion  of  this 
program.  The  Raytheon  effort  emphasized  the  evaluat  ion  two-dimensional 
adaptive  arrays  ir.  the  real  HF  environment.  The  primary  task  was  to  modify 
an  existing  HF  phased  array  facility  in  Hudson,  Colorado  so  that  data  could  be 
'ie.cted  on  a large  two-dimensional  array  in  a format  suitable  for  use  by  Ray- 
theon and  by  others  for  adaptive  proces  the  principal  modifications 

included  the  construcdon  of  a 3'  n*  cross  array  of  vertical  monopoles 

and  the  addition  of  quadrat  outputs  to  the  existing  32-channel  receiver  and 
data  acquisition  systet  . Additional  Raytheon  tasks  included  the  development 
of  software  (for  controlling  the  data  collection  system,  for  processing  and  dis- 
playing the-  field  data,  for  producing  simulated  field  tapes  for  specific  array 
geometries  and  model  environments),  data  collection  under  a variety  of  condi- 
tions, and  analysis  (numerical  simulations  to  guide  the  design  of  the  antenna 
array  and  data  reduction  to  evaluate  the  performance  improvement  obtained 
through  adaptive  processing). 

The  principal  result  of  this  work  is  the  verification  that  the  HF  environ- 
ment will  in  fact  permit  most  of  the  performance  improvement  anticipated 
from  idealized  numerical  simulations  of  adaptive  processing  applied  to  large 
two-dimensional  arrays.  Noise  f'.oor  reductions  of  order  30  dB  (compared 
to  conventional  processing)  were  realized  over  most  of  the  radio  sky  under 
conditions  in  which  a single  strong  interferer  was  present  (10  dB  element  sig- 
nal to  noise  ratio).  Viewed  as  a technique  for  producing  improved  two- 
dimensional  maps  of  the  radio  environment  (e.  g.  , radio  location  or  noise 


survey  applications)  rather  than  for  improving  signal  to  noise  ratio,  the  maps 
produced  by  adaptive  beamfo r rning  exhibit  superior  angular  definition  and 
appa  ri  n t freedom  f-om  sid.-lobe  responses.  That  is,  tin  adaptive  processor 
,s  in  deconvolving  the  antenna  pattern  from  the  m«  asurements. 

T,"‘  MF  ‘ nviionment  is  found  to  change  rapidly  enough  so  that  if  a one 


second  block  of  data  is  processed  with  adaptive  weights  computed  for  the  pre- 
ceding one  s<  corn!  data  block,  approximately  10  dR  of  the  noise  floor  improve- 
mc  nt  attainable  by  uaing  opfmum  weights  is  lost.  It  appears  that  the  direct 

' inversion  method  for  implementing  adaptive  processing  will  in  practice 
require  a new  inversion  for  each  data  block  but  it  has  not  yet  been  determined 
bow  long  these  data  blocks  ran  be  made  without  incurring  an  unacceptable  sig- 
nal to  noise  ratio  penalty. 


It  was  found  that  the  output  of  the  adaptive  processor  was  typically  20  dB 
less  than  that  of  the  conve  ntional  processor  when  both  were  pointed  in  the  source 
rect  ion,  an  unexpected  effect  (ideally  they  would  be  equal)  attributed  to  the 
departure  of  the  wavefront  of  the  incident  signal  from  ar.  ideal  plane  wave. 

Hus  departure  is  believed  due  to  ionospheric  multipath  but  terrain  scattering 
and/or  residual  calibration  errors  in  the  receiving  system  cannot  be  ruled  out. 

The  result  of  the  array  geometry  study  is  to  suggest  low-redundancy  sparse 
arrays  as  promising  candidates  for  adaptive  antenna  applications  in  which  both 


angular  resolution  and  signal  to  noise  ra'io  are  considerations. 

It  is  recommended  that  the  processing  of  the  data  base  already  on  tape  be 
completed  so  as  to  explore  the  impact  of  variable  data  block  length  on  signal  to 
noise  ratio,  to  explore  the  consequences  cf  simultaneous  adaptivity  in  azimuth, 
elevation,  and  frequency,  and  to  explore  the  performance  of  the  adaptive 


processor  in  complex  interference  environments. 
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1.0  INTRODUCTION 


U- 


1. 1 object:  e s 

The  DARPA  HF  Adaptive  Antenna  program  has  been  a three  component 
program  involving  the  Stanford  Research  Institute  (SRI),  the  University  of 
Colorado  (U  of  C)  and  Raytheon.  The  overall  objective  of  the  program  has 
>ecn  to  determine  the  applicability  of  known  adaptive  processing  techniques 
to  Hi  antenna  arrays.  The  motivation  for  the  program  has  come  from  the 
impressive  improvements  in  array  performance  obtained  ly  applying  adaptive 
processing  to  sonar,  seismic,  and  microwave  systems.  The  need  for  an  HF 
program  arises  because  the  benefits  derived  from  adaptive  processing  depend 
strongly  upon  the  characteristics  of  the  noise  and  interference  environment 
(e.  g.  , temporal  ant!  spatial  stability):  consequently,  performance  improve 

ments  obtained  at  microwave  frequencies  arc  not  necessarily  good  indicators 
of  what  can  be  expected  at  HF  frequencies.  To  obtain  the  necessary  data  base 
at  HF  and  to  test  the  applicability  of  adaptive  algorithms  developed  with  other 
environments  in  mind,  the  DARPA  program  was  structured  in  the  following 
' way:  SRI  modified  the  receiver  system  for  the  WARF  antenna  so  as  to  pro- 

vide data  for  adaptive  processing  from  a large  one -dimensional  array, 
Raytheon  made  modifications  to  the  DARPA  HF  phased  array  facility  at  Hud- 
son, Colorado  so  as  to  provide  data  for  adaptive  processing  on  a large  two- 
dimensional  array,  and  U of  C provided  algorithm  research  and  processing 
expertise.  An  interim  report  on  the  Raytheon  component  of  the  program  was 
published  in  March  1974:  (Reference  1)  the  material  preserved  here  is  the 

final  report  on  this  effort. 
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1 .2 


SCOPE  OF  THE  EFFORT 


The  primary  task  (70%  of  the  effort)  was  to  modif/  the  Hudson  phased 
array  facility  so  that  data  could  be  enllected  on  a large  two-dimensional  array 
in  a format  suitable  for  use  by  ourselves  and  by  others  for  adaptive  process- 
ing. Two  principal  modifications  were  required:  First,  the  receiver  system 

was  modified  to  provide  simultaneous  quadrature  samples  for  each  element  of 
a 52  element  array.  The  original  receiver  system  consisted  of  a bank  of  32 
receiver  modules,  each  with  a single  real  output  and  a single  sample  and  hold 
ircuit.  Quadrature  outputs  were  provided  by  fabricating  an  identical  set  of 
12  receiver  modules,  shifting  the  phase  of  the  local  oscillator  (HFO)  for  this 
set  by  90°  from  the  first  set,  and  using  one  receiver  module  from  each  set  'or 
each  antenna.  This  provided  the  complex  samples  required  in  adaptive  proc- 
essing.  The  second  facility  modification  was  to  build  a cross  array  of  32 
vertical  monopole  elements  with  base-  amplifiers.  The  original  facility  array 
was  a circular  array  with  horizontal  crossed  dipole  elements  intended  for  diag- 
nosing the  overhead  ionosphere.  This  array  was  not  suitable  for  studying  the 

i 

low  elevation  angle  signals  of  primary  interest  in  adaptive  antenna  applications. 

The  second  major  task  (15%  of  the  effort)  was  the  development  of  soft- 


ware for  three  purposes: 

I 

1.  field  software  for  controlling  the  data  acquisition  system 
and  for  providing  on-site  diagnostic  information. 

2.  analysis  software  for  processing  and  displaying  the  eld  data 

using  conventional  and  adaptive  two-dimensional  beamforming 

I 

algorithms.  / 

3.  simulation  software  to  generate  data  tapes  for  model  signal 
environments/ and  for  model  array  configurations  to  use  in 


/ 


testing  thr  analysis  software  and  for  studying  the  effects  oi  array 
configuration  on  processor  performance. 

The  third  task  was  data  collection  (5%  of  the  effort)  which  involved  pre- 
paring and  manning  the  Hudson  phased  array  site  for  three  flays  01  data  taking 
( l '.vo  days  in  November,  1973  and  one  day  in  January,  1974).  Data  were  collected 
with  SKI  providing  cw  transmissions  at  several  different  power  levels  anti 
frequencies.  Data  were  also  collected  in  a variety  of  uncooperative  interfer- 
ence environments  such  as  the  amateur  radio  bands  and  the  international 
broadcast  bands. 

The  fourth  task  was  data  analysis  (lO^o  of  the  effort).  The  analysis 
involved  the  reduction  of  selected  portions  of  the  data  base  with  conventional 
and  adaptive  proccssi  ng  and  the  evaluation  of  the  results  in  determine  such 
things  as  the  stability  of  the  environment  and  the  noise  floor  reduction 
achieved.  This  task  also  included  conjidcring  the  effect  of  array  geometry 
on  adaptive  processar  performance. 

1.  3 ORGANIZATION  OF  THE  REPORT 

This  report  is  written  for  the  reader  who  is  familiar  with  HF  systems 
and  the  problems  they  encounter  when  operated  in  the  real  world:  a back- 
ground in  adaptive  processing  is  not  assumed.  The  presentation  begins  (Sec- 
tion 2.0)  with  a brief  statement  of  the  theory  and  the  processing  algorithms  on 
which  our  work  has  been  based.  A series  of  computer  simulations  is  then 
presented  (Section  3.0)  to  illustrate  the  operation  of  adaptive  processors 
under  ideal  conditions,  to  document  the  research  which  has  been  done  on  the 
properties  ol  low-redundancy  sparse  arrays  when  used  with  adaptive  beam- 
formers,  and  to  provide  the  basis  for  comparing  the  theory  with  the  field 
measurements.  Section  4.0  includes  a description  of  the  two-dimensional 


3 


antenna  array  used  in  this  program  and  a discussion  of  the  data  base  collected. 
An  analysis  of  selected  portions  of  this  data  base  is  given  in  Section  5.0.  The 
results  of  the  Raytheon  component  of  the  program  are  summarized  in 
Sectir n 6.  0. 
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2.  0 THEORETICAL  BACKGROUND 


The  purpose  of  this  section  is  to  give  a brief  statement  of  the 
mathematics  which  has  been  used  in  the  software  codes  for  conventional  and 
adaptive  processing.  The  mathematics  is  particularly  simple  and  compact  in 
matrix  notation  and  this  form  of  the  equations  will  be  used  for  convenience. 

T his  section  also  includes  a discussion  of  the  effect  of  receiver  errors  on 
i leptive  beamformers,  in  preparation  for  the  analysis  of  experimental  data 
to  be  presented  in  Section  'i.  0. 

2.  1 BE/MFORMER  OPERATION 

Bcamforming  is  the  process  o(  estimating  the  sigi.al  power  arriving 
from  a given  direction  and  amounts  to  summing  the  signals  from  each  array 
element  after  applying  an  appropriate  set  of  complex  weights.  The  process 
is  illustrated  schematically  in  Figure  1 and  is  the  same  operation  whether 
conventional  or  adaptive  processing  has  been  used  to  determine  the  weights. 

It  is  assumed  that  the  array  has  M elements  and  that  the  receiver  output  from 
each  clement  is  sampled  simultaneously  to  give  a set  of  M complex  samples, 
x t x xJVf  These  samples  arc  weighted  with  a set  of  complex  weights 

to  form  a new  data  set  (x^,  XMWM)  and  this  nCW  data  SCt  iS 

summed  to  give  a complex  number,  y,  as  the  bcamformer  output.  If  the  set 
of  weights  is  written  as  a column  vector  W and  the  set  of  data  samples  as  a 
column  vector  X,  then  the  beamformcr  operation  can  be  written  in  matrix 
notation  as 

y ^ W':X  0) 
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whc  re  W - is  the  conjugate  transpose  of  W (i.  e.  , a row  vector  with  elements 
equal  to  the  complex  conjugate  of  the  elements  in  W).  This  gives  the  beam- 
form,  r output  for  one  time  sample  of  the  array  outputs.  In  this  report  this 

, Ss  is  always  ••epetei  lor  (>‘\  time  samples  of  the  array,  the  power  out- 
put compute, I for  earn  sample,  ami  the  results  averaged  to  give  the  average 
now,  r arriving  from  the  direction  of  interest. 


I . i CONVENTIONAL  WEIGHTS 

ta  conventional  beamforming,  the  computation  of  the  weights  to  form  a 
beam  in  a gi\en  direction  begins  by  computing  the  signal  expected  at  each 
element  of  the  array  for  a plane  wave  of  unit  amplitude  arriving  from  the 

desired  direction.  Tnis  set  of  complex  numbers  (vj,  v£ vm*  is  writton 

as  a column  v,  ctor  V and  this  is  called  the  steering  vector  for  that  direction 
In  conventional  beamforming  the  weight  vector  is  set  equal  to  the  steering  vec  - 

tor  a nd  cons,  ,|iient ly 

W*  V'v 

for  use  in  the  beamformer  of  equation  (1).  A constant  of  proportionality  could 
be  added  to  the  righ.  side  of  this  equation  without  affecting  the  antenna  pattern 
or  the  SNR  but  we  elect  to  choose  this  constant  to  be  unity  in  order  to  keep  the 
expression  as  simple  as  possible.  With  this  choice,  the  output  of  the  conven- 
tional beamformer  will  be  M units  when  steered  to  look  in  the  direction  ot  a 

plane  wave  of  unit  amplitude  at  each  element. 

If  the  background  noise  environment  is  omni-directional  (a  "white  1 
angular  noise  spcctium)  the  conventional  weights  are  known  to  give  the  maxi- 
mum signal  to  noise  ratio  (SNR)  out  of  the  beamformer.  In  practice  this 
condition  is  seldom  if  ever  met  and  in  this  case  (a  "colored"  angular  noise 
spectrum)  the  optimum  weights  must  be  adapted  to  match  the  environment  if 
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maximum  SNR  is  tc  be  obtained. 


2.  3 ADAPTIVE  WEIGHTS 

Tin  first  step  in  computing  the  adaptive  weights  is  to  compute  the 
average  covarianee  matrix  for  the  array.  This  covariance  matrix,  F,  is  an 
M by  M matrix  made  up  of  the  average  covariance  between  all  pairs  of  elements 
(<  . g.  , fj  j is  the  covariance  between  the  first  element  and  the  seventh).  This 
matrix  characterizes  the  environment  ove.  the  period  of  interest.  The  conjug- 
al* transpose  of  the  weight  vector  is  then  given  (Griffiths,  Reference  2)  by 

W = b V " F"1  (3) 

wlu  re  b is  a constant  which  normalizes  the  set  of  weights  in  a convenient 
manner,  V is  the  conjugate  transpose  of  the  steering  vector,  and  F ^ is  the 
inverse  of  the  covariance  matrix.  It  is  convenient  to  choose  the  constant  b so 
that  the  output  of  the  adaptive  bcamformer  will  equal  the  output  of  the  conven- 
tional beamformer  when  both  are  steered  in  the  direction  of  a plane  wave  of 
unit  amplitude  at  each  element.  The  appropriate  constant  (Capon,  Reference 
3)  is 

b = VI/ V*  F’ 1 V (4) 

Equations  (1),  (3),  and  (4)  constitute  the  adaptive  beamformer  as  implemented 
in  this  report. 

2.  4 RECEIVER  ERRORS 

In  any  practical  beamformer  the  received  data  samples  are  somewhat 
in  error  because  the  M data  channels  have  residual  gain  and  phase  shift  differ- 
ences between  chanrels.  The  effect  of  these  errors  has  been  discussed 
extensively  for  conventional  beamforming.  The  effect  of  errors  in  the  data 
samples  or  in  the  weighting  function  of  adaptive  processors  has  also  received 
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some  attention.  McDonough  (Reference  4)  has  provided  some  insight  into 
effects  of  such  errors,  and  relates  his  results  to  those  of  other  authors.  He 
points  out  that  in  conventional  beamforming,  the  principal  effect  of  errors 
is  in  the  "sidelobe1  region  of  an  antenna  pattern.  With  adaptive  or  optimum 
processors,  a degrading  effect  on  the  resolution  properties  near  a discrete 
source  is  noted. 

As  McDonough  notes,  the  effect  of  errors  in  the  receiving  channels 
will  vary  depending  on  how  the  errors  appear  in  the  adaptive  beamformer. 

1 lms,  it  is  desirable  to  consider  the  effect  of  errors  in  the  manner  in  which 
tiny  are  most  likely  introduced  in  the  system  being  evaluated. 

For  the  experimental  data  discussed  in  this  report,  the  environment 
is  essentially  a discrete  interfering  source  and  isotropic  noise.  With  channel 
to  channel  gain  ar.d  phase  unbalance,  the  residual  errors  cause  a modification 
of  1 he  oils  • rved  covariance  matrix  relative  to  the  true  covariance  matrix  with 
no  errors.  I’he  per  urbed  covariance*  matrix  is  the  one  used  to  form  the  adap- 
tive weights.  The  r -spouse*  of  the*  processor  will  be  el?  ferent  from  the  response 
with  no  errors.  The  effect  of  these  errors  on  the  adaptive  processor  is 
analyzed  below:  simulations  are  presented  to  check  the  analysis  for  the 
representative  case  of  a discrete  source  in  isotropic  noise. 

<£.4.1  Error  Analysis 

The  samples  trom  the  environment  in  the  absence  of  errors  consist  of 
a noise  power  c f X watts  at  each  element  (uncorrelated  from  element  to  ele- 
ment) and  a plane  wave  signal  of  S watts  at  each  element.  The  covariance 
matrix  without  errors  is 

F = NI  + SVV*  (5) 
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here  1 is  the  unit  diagonal  matrix,  and  V is  the  steering  or  signal  vector  tor 
the  discrete  source,  corresponding  to  the  direction  of  incidence  of  the  plane 
wave.  The  principal  effect  of  channel  to  channel  errors  is  to  produce  a cor- 
rupird  signal  vector,  whose  elements  are 

small  random  gain  and  phase  changes,  c.. 


\V 


nominally  those  of  V,  but  with 


i.  '!  ± E/t)  e 


±E/l 


(6) 


re  E is  a small  number  reflecting  th-  error  in  amplitude  or  in  phase. 

1 h.  present  analysis  and  simulation  assumes  a pseudorandom  binary  error 
model  with  independent  phase  and  amplitude  errors.  It  would  be  straight- 
forward to  extend  the  model  to  Gaussian  errors  with  different  magnitudes  of 


am 


plitude  and  phase  as  given  by  McDonough  (Reference  4).  The  present 
analysis  continues  by  applying  this  error  to  the  discrete  source.  Further 
analysis  and  simulation  have  shown  that  the  effect  of  cor. option  of  isotropic 
noise  by  these  errors  is  negligible  in  comparison  to  that  of  ll.e  discrete  source, 
when  the  source  power  S is  equal  to  or  greater  than  the  noise  power.  N. 

The  observed  covariance  matrix  which  contains  the  effect  of  these 


errors  is 


(7) 


Fc  = NI  + svcVc 

When  the  adaptive  beamformcr  is  steered  to  the  direction  of  the  inter- 


ference,  the  reciprocal  of  the  output  power 


is 


p _l  = MN 
v 


1 + (S/N)  V*  Vc  -(S/MN)  \V^[ 


1 + (S/N)  Vn  V 


e e 


(8) 


This  expression  results  in  an  output  power  of  M2S  watts  duo  to  the  plane  wave 
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input  where  E 0.  Next  it  is  noted  th3.t  when 

(S/N)  V V -(S/MN)  jVe  V j >>  1 (9) 

tin-  expression  (H)  for  the  output  power  Py  is  independent  of  the  input  sienal 
powor  S.  This  means  there  is  a "saturation"  effect  in  the  output  of  an  adap- 
tive' beamformer  when  the  conventional  output  SNR,  equal  to  SM  N,  is  large 
i nough. 

I'o  evaluate  the  criterion  for  saturation,  the  pseudorandom  error 

mo  lei  is  used  to  compute  the  values  of  ani  |v^  v]  . In  estimating 

these  quantities,  the  approximation  will  provide  representative  values  for  Py 

and  the  criterion  for  saturation.  Using  the  values  for  V"  M (1  + F7  / A)  and 

2 ..2,.  -2 


V V 2w  m2(1  -t  / 4),  (d)  becomes 

V 


E2/2>>  (SM/N) 


(10) 


and  (H)  becomes 


Thus  the  criterion  for  saturation  equation  (10),  is  tnat  the  error  magni- 
tudes must  be  greater  than  the  reciprocal  square  root  of  the  conventional  output 
signal  to  noise  power  ratio.  When  this  criterion  is  met,  the  output  power,  (11) 
is  a constant  regardless  of  any  increase  in  signal  power  S. 

This  analyses  can  he  extended  to  consider  look  directions  away  from 
the  interfering  source.  In  those  cases,  the  output  power  is  not  significantly 
affected  as  long  as  the  error  is  small  compared  with  unity.  Thus  the  adap- 
tive beamformer  is  less  sensitive  to  channel-to-channel  errors  than  the  con- 
ventional beamformer  in  the  "sidelobe"  regions. 
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2 . *1 . 2 Simulation 


Simulations  wort  performed  using  a pseudorandom  error  code,  ont 
for  magnitude  and  another  for  phase,  with  a 32  element  limar  array.  I igur< 

2 shows  that  for  a fixed  element  SNR,  the  amount  by  which  the  beamformer 
unci i rest imates  the  signal  level  is  small  tor  small  receiver  errors,  in*  rcases 
slowly  with  increasing  receiver  error  at  first  and  then  a break-point  is 
reached  after  which  the  underestimate  increases  much  more  rapidly  with 
receiver  error.  As  the  error  increases,  the  break-point  is  reached  sooner 
tor  high  element  SNR  than  for  low,  and  consequently  an  adaptive  system  with 
.a  given  set  of  receiver  errors  will  give  better  estimates  of  the  signal  power 
for  weak  signals  than  for  strong.  Once  the  break-point  has  been  passed,  the 
beamformer  output  saturates,  in  that  further  increases  in  element  SNR  pro- 
duce no  further  increase  in  the  output.  For  example,  for  an  RMS  gain  error 
of  10  percent,  a tO  dB  increase  in  element  SNR  from  20  d B to  30  dB  results 
in  an  additional  10  dB  underestimate  of  the  signal  level,  just  compensating  for 
the  increase  in  element  SNR  and  causing  the  output  to  remain  constant. 

It  has  been  assumed  in  this  discussion  that  the  errors  are  caused  by 
receiver  gain  and  phase  unbalance,  causing  an  incident  plane  wave  to  look 
irregular  to  the  processor.  Alternately,  the  errors  can  be  thought  of  as 
being  introduced  by  irregularity  in  the  ionosphere,  by  irregularity  in  the 
ground  from  clement  to  element,  or  by  irregularity  caused  by  scattering  from 
local  structures.  Whatever  the  cam.  , the  adaptive  processor  will  underestimate 
the  signal  power  if  the  incident  wave  front  does  not  appear  plane  to  the  processor: 
the  curves  of  Figure  2 allow  one  to  estimate  how  serious  this  effect  will  be. 
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3.0  SIMULATION  RESULTS 


The  algorithms  discussed  in  Section  2.0  for  processing  two-dimensional 
array  data  conventionally  and  adaptively  were  implemented  in  software,  llong 
with  graphical  routines  for  displaying  the  results.  This  software  package  vas 
wrilt«  ti  to  accept  data  tapes  from  the  field  and  has  been  usee!  to  reduce  the 
field  data  discussed  in  Section  5.0.  A separate  software  package  was  also 
written  to  generate  simulated  data  tapes  for  model  signal  environments  and 
model  array  geometries.  That  is,  the  intensity,  azimuth,  elevation,  and 
Doppler  shift  of  a finite  number  of  point  sources  could  be  specified  along  with 
Hie  intensity  of  an  isotropic  background  noise  level  and  the  number  and  place- 
ment of  antenna  elements  in  a one-  or  two-d'mensional  array.  The  simulated 
data  tape  contained  the  time  sampled  outputs  from  each  element  of  the  array 
which  would  be  expected  under  these  idealized  conditions.  These  data  tapes 
could  then  be  processed  as  though  they  were  field  tapes.  The  simulations  were 
used,  1)  to  verify  and  debug  the  processing  software  (by  simulating  special 
cases  in  which  the  rerrect  processor  output  was  known),  2)  to  study  the  effect 
of  array  geometry  on  processor  performance  (to  decide  how  to  make  the  best 
use  of  the  32  elements  available  for  the  collection  of  field  data)  and  3),  to  pro- 
vide test  cases  with  which  to  judge  whether  or  not  the  properties  of  adaptive 
processors  as  predicted  from  idealized  models  were  realized  in  practice  with 
’•cal  field  data.  A series  of  test  cases  under  idealized  conditions  are  given 
below  (Section  3.  11  to  illustrate  the  properties  expected  of  adaptive  processors. 
The  results  of  our  array  geometry  study  are  given  in  Section  3.2. 


Preceding  page  blank 
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P KOPERTILS  of  ideal  adaptive  processors 


3.  1 


3.1.1  Adaptive  Array  Antenna  Patterns 

This  section  illustrates  the  interference  rejection  principle  of  adaptive 
Ik  miforiners  by  examining  the  way  in  which  the  antenna  pattern  changes  from 
one  look  direction  to  the  next.  Stated  briefly,  an  adaptive  heamformer  affects 
in  improvement  is  SP'R  (as  compared  to  a conventional  heamformer)  by  adap- 
tively tailoring  the  antenna  array  weights  (both  in  amplitude  and  phase)  in 
< rdt  r to  set  nulls  in  the  corresponding  array  pattern  in  the  direction  of  dis- 
crete interferers  A suitable  array  gain  is  simultaneously  maintained  in  the 
desired  look  direction  so  that  SNR  is  maximized. 

Consider  the  seven  element  linear  array  shown  in  Figure  3.  The  ele- 
ments are  isotropic  and  spaced  by  half  a wavelength.  The  simulated  signal 
environment  consists  of  three  discrete  sources  situated  at  -40,  1?,  and  30 
degrees  off  boresight  with  relative  strengths  of  -20,  0,  and  0 dB,  respectively. 

1 he  strength  of  the  omni-directional  noise  is  -30  dB  (all  signal  strengths  meas- 
ured at  an  individual  element  before  any  array  or  processing  gain).  Working 
against  this  environment,  the  adaptive  weights  have  been  computed  for  look 
directions  of  -4  and  -40  degrees:  the  magnitudes  of  the  complex  weights  for 
these  two  look  directions  are  shown  as  a function  of  element  number  at  the 
top  of  Figure  4.  The  c ir responding  patterns  of  the  array  for  these  two  look 
directions  are  shown  at  the  bottom  of  this  figure.  By  the  "pattern"  for  a given 
look  direction  we  mean  a plot  of  the  power  output  of  the  heamformer  (with  the 
adaptive  weights  for  that  look  direction  set  in  and  held  fixed)  versus  the  angle 
to  a point  source  of  unit  strength.  From  patterns  such  as  these  the  following 
observations  can  be  drawn: 

1.  The  array  pattern  nulls  track  the  directions  of  the  discrete 
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igure  3.  Model  Environment  Used  for  Simulations  and  a Seven 

Element,  Uniformly  Spaced  Array  of  Isotropic  Radiators 


interferers  as  the  desired  look  direction  is  changed. 


L.  The  depth  of  the  antenna  pattern  nulls  are  related  to 

(he  signal  strengths  of  the  corresponding  discrete  inter- 
ferer s . Deeper  nulls  are  placed  on  the  stronger  interferers. 

i.  Hu  adaptive  array  patterns  typically  have  a response  maximum 
near  the  desired  look  direction. 

1 I he  . daptive  array  response  patterns  are  not  "super- 
directive".  Their  "main  beam"  widths  are  similar  to 
those  of  a uniform  filled  conventional  array.  Also,  the 
very  high  (e.  g.  , 10^)  weight  magnitudes  characteristic 
of  super  gain  schemes  have  not  been  observed. 

S.  The  adapted  weight  values  increase  for  look  directions  in 
the  vicinity  of  discrete  signal  components  of  high  SNR. 

3.1.2  Adaptive  '’e-sus  Conventional  X^aps  of  the  Environment 

The  concept  of  an  array  pattern,  while  still  valid,  is  of  limited  useful- 
ness in  describing  the  performance  of  adaptive  processors  used  to  scan  the 
environment  (say  io  map  the  location  of  interferers  or  to  search  for  weak 
target  echoes).  In  contrast  to  the  patterns  of  conventional  beamformers 
which  remain  constant  (in  sine  space)  as  the  look  direction  is  scanned,  the 
pattern  of  an  adaptive  arr?y  is  different  for  each  look  direction,  there  may 
or  may  not  be  a "main  beam"  and  it  may  or  may  not  point  in  the  desired  look 
direction,  and  the  peak  sidelobe  level  of  the  pattern  for  ar.y  given  look  direction 
is  a poor  indicator  of  how  serious  the  effects  of  strong  interferers  will  be  when 
the  beam  is  steered  elsewhere.  It  will  be  more  useful  to  think  of  the  adaptive 
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processor  as  an  estimator  of  the  power  arriving  from  a given  direction  and  to 
ignore  the  details  oi  the  pattern  formed  in  order  to  obtain  that  estimate.  Ihis 
sc*,  lion  will  compare  the  quality  of  that  estimate  for  conventional  and  adaptive 
beamformers. 

Kieur.  S compares  the  sky  maps  (signal  intensity  versus  angle)  ob- 
tain' d by  si  .inning  t lie  environment  with  a conventional  beamformer  (top;  and 
an  adaptive  beamformer  (bottom).  The  array  and  the  environment  are  the 
•ame  as  illustrate*!  in  Figure  3.  The  conventional  beamforming  was  done 
Kb  nutorm  weights  tno  taper).  Arrows  indicate  the  position  of  the  dis- 
rnti  signal  sources.  The  conventional  map  suffers  from  two  problems:  the 
two  strong  sources  are  not  well  resolved  and  the  existence  of  the  weak  signal 
is  masked  by  side  lobe  responses  from  the  two  strong  signals.  These  prob- 
1,  ms  are  not  present  on  the  adaptive  map.  The  two  strong  sources  are  well 
resolved  and  the  weak  source  is  clearly  visable.  The  adaptive  processor  has 
miU|,  )><  tter  use  of  the  information  available  from  the  array  elements  in  estim- 
ating the  characteristics  of  the  environment.  The  improvement  can  be 
described  by  saying  that  the  effective  bearmvidth  of  the  array  has  been  nar- 
rowed and  that  the  affective  sidelobe  level  of  the  array  has  been  reduced,  if 
it  is  kept  in  mind  that  the  array  pattern  for  any  given  look  direction  does  not 

exhibit  these  proper*ies  (e.  g.  , Figure  4). 

A second  comparison  between  sky  maps  produced  with  conventional 
and  adaptive  processors  is  shown  in  Figure  6.  In  this  case  a 40  dB  Dolph 
taper  has  been  applied  to  the  conventional  array  in  order  to  suppress  side- 
lobes.  It  is  quite  effective  in  doing  so  but  at  the  cost  of  broadening  the  main 
beam  and  thus  making  the  problem  of  resolving  the  two  strong  sources  worse. 
The  adaptive  map  is  shown  as  before  for  comparison.  Not  only  are  the  sources 
better  resolved  in  <he  adaptive  display  but  the  noise  floor  between  sources  is 
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Figure  S.  Comparison  between  the  Performance  of  a Conventional 
Beamformer  with  Uniform  Weights  (top)  and  an  Adaptive 
Beamformer  (bottom). 
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significantly  lower. 

The  improvement  in  SNR  achieved  by  applying  adaptive  processing 
instead  of  conventional  processing  (with  uniform  weights'  to  the  array  and 
environment  shown  in  Figure  3 is  plotted  as  a function  of  look  direction  in 
Figurt  7.  This  is  the  difference  between  the  two  curves  shown  in  Figure  5 
and  illustrates  that  the  improvement  obtained  varies  great!'  from  one  look 
direction  to  another.  No  improvement  is  obtained  when  looking  directly  at 
disc  rele  sources  r.nd  little  improvement  is  obtained  in  directions  such  that 
llu  nulls  between  iho  sidelobes  of  the  conventional  pattern  straddle  the  strong 
interferers . 

These  examples  illustrate  the  features  which  characterize  adaptive 
processors  and  make  them  potentially  attractive  for  HF  rar.ar  and  communi- 
cations applications.  These  examples  are  of  course  special  cases  and  care 
must  be  exercised  when  extrapolating  to  other  situations  because  the  perform- 
ance achieved  will  depend  upon  both  the  array  and  the  environment. 

3.  i ADAPTIVE  PROCESSORS  WITH  SPARSE  ARRAYS 

Simulations  have  been  used  to  study  the  impact  of  element  placement 
(array  geometry)  on  the  performance  of  an  adaptive  processor.  The  primary 
motivation  for  this  study  was  the  need  to  choose  a specific  array  geometry  for 
installation  in  the  field  to  collect  the  program  data  base.  The  objectives  of 
the  data  collection  were 

1.  to  map  the  environment  so  as  to  define  the  strength  and 
position  (azimuth  \nd  elevation)  of  interferers, 

and 

2.  to  evaluate  the  improvement  (e.g.,  noise  flooi  reduction) 
achieved,  as  a function  of  position  relative  to  the  interferers, 
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Figure  7.  Impiovement  in  SNR  through  Adaptive  Processing  Rather 
than  Conventional  Processing  with  Uniform  Weights. 


by  processing  the  array  adaptively  rather  than  conventionally. 

The  constraint  was  that  only  32  elements  (and  their  associated  receivers  and 
pI.(>  cssing  equip.m  n»)  . ould  be  afforded  within  the  program  budget  and  the 
problem  was  to  distrdmte  these  in  two-dimensions  so  as  to  optimize  tin  pe  r - 
formant  e of  the  array.  The  choice  then  was  between  : mall  filled  array  or 

a large  sparse  array.  In  conventional  beamforming  the  small  array  gives 
better  sidelobe  performance  but  poorer  angular  resolution. 

It  was  first  recognized  by  radio  astronomers  that  both  high  angular 
resolution  and  low  sidelobe  levels  could  be  achieved  at  the  same  time  for  a 
special  class  of  sparse  arrays  (low  redundancy  arrays)  if  special  processing 
techniques  were  employed  (Reference  4).  The  basis  for  mis  approach  is  the 
fact  that  the  intensity  of  the  angular  spectrum  (the  sky  map)  can  be  otained 
either  by  taking  the  Fourier  transform  of  the  aperture  distribution  and 
squaring  the  magnitude  of  the  result  (this  is  the  software  implementation  of 
conventional  processing)  or  by  first  computing  the  autocorrelation  function 
across  the  aperture  (correlation  coefficient  versus  interelement  spacing)  and 
Fourier  transforming  this.  The  significance  of  this  is  that  it  shows  that  all 
the  sky  map  information  of  interest  is  contained  in  the  measured  autocorrela- 
tion function  for  the  environment:  each  point  on  the  autocorrelation  function 

represents  a different  inter-element  spacing  in  the  array  and  (for  averaging 
times  long  enough  so  that  the  incident  signals  become  incoherent  with  each 
other)  the  result  for  a given  spacing  is  the  same  no  matter  which  pair  of 
elements  at  that  spacing  is  used.  This  means  that  elements  of  a filled  array 
can  be  removed  so  long  as  all  the  spacings  present  in  the  original  filled 
array  still  occur  at  least  once.  The  elements  saved  in  this  way  can  be  added 
outside  the  boundary  of  the  original  filled  array  so  as  to  make  available  large 
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int<  r-element  spacings  not  present  at  all  in  the  original  array.  Ideally, 
array  configurations  could  be  found  for  which  all  spacings  occur  only  onc< 
such  arrays  are  called  non- redundant  or  zero- redundancy  arrays  but  exist 
only  for  arrays  of  four  elements  or  less.  For  larger  numbers  of  elements 
some  spacings  occur  more  than  once  and  these  are  called  low-redundancy 
arrays.  The  configuration  which  gives  the  lowest  number  of  redundant  spac- 
ings for  a given  number  of  elements  is  called  a minimum  redundancy  array. 

In  my  event,  the  processing  procedure  is  to  compute  the  autocorrelation 
1 unction  for  the  array,  weight  this  function  to  control  sidelobes  (as  one  would 
weight  the  physical  aperture  in  conventional  processing  to  control  sidelobes), 
and  Fourier  transform  the  weighted  autocorrelation  func*ion  to  obtain  the  sky 
map.  The  key  point  is  that  it  is  the  autocorrelation  function  for  the  array 
which  must  be  filled,  not  the  physical  aperture:  the  physical  array  can  be 
made  sparse  while  maintaining  low  sidelobe  levels  by  aop’ying  weights  to 
the  autocorrelation  function  rather  than  to  the  physical  aperture.  The  array 
can  be  thinned  until  further  element  removal  would  cause  a "hole1  in  the  auto- 
correlation function  (a  spacing  not  represented  in  the  array). 

If  mapping  of  the  interference  environment  had  been  the  sole  purpose 
of  this  program,  a sparse,  low-redundancy  array  would  have  been  the  obvious 
choice  with  processing  and  sidelobe  control  done  non-ada ptively  in  the  auto- 
orrelation  domair.  It  was  not  obvious  that  such  arrays  would  also  be  suitable 
for  adaptive  processing  but  it  seemed  plausible  that  this  would  prove  to  be  the 
case  for  the  following  reason:  adaptive  processors  have  the  property  that  they 
try  to  adjust  the  pattern  of  the  array  so  that  the  sidelobes  don't  fall  on  strong 
interferers.  Even  when  used  with  filled  arrays  capable  of  low  pattern  side- 
lobe  levels,  adaptive  processors  characteristically  produce  patterns  with 
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high  sidelobes  because  it  is  often  preferable  from  an  SNR  viewpoint  to  have 
high  sidelobes  pointing  in  quiet  directions  in  order  to  obtain  deep  nulls  in  just 
the  right  direction  to  reject  strong  discrete  interferers.  That  is,  pattern 
flexibility  is  more  important  than  peak  sidelobe  level.  The  fact  that  sparse 
d r rays  have  high  oeak  sidelobes  does  not,  therefore,  necessarily  disqualify 
them  as  effective  adaptive  arrays,  particularly  when  the  elements  saved  by 
making  the  array  sparse  are  used  to  make  the  overall  dimensions  of  th  array 
larger  and  thus  to  give  the  adaptive  processor  new  and  longer  baselines  to 

| 

work  with. 

A series  of  simulations  were  done  to  study  the  performance  of  low 
redundancy  spars"  arrays  under  ac  aptive  processing.  The  first  case  illus- 
* trated  compares  the  performance  of  a filled  seven  element  array  to  a sparse 

i 

four  element  array  cf  the  same  aperture.  The  filled  array  and  the  environ- 
ment is  the  same  as  that  shown  in  Figure  3.  The  sparse  array  is  formed  by 
removing  elements  number  3,  5,  and  6:  the  resulting  array  has  all  the  inter- 
element spacings  of  the  filled  array.  Figure  8 compares  the  performance  of 
this  array  (top)  with  the  performance  of  the  fi1  led  array  (bottom)  under  adap- 
tive processing.  The  filled  array  does  a better  job  of  maintaining  a uniformly 
low  noise  fuor  but  the  source  resolution  achieved  with  both  arrays  is  compar- 
able (not  surprising,  since  the  ape^ure  is  the  same)  and  the  sparse  array 
with  adaptive  processing  is  much  superior  to  a filled  array  of  the  same  aper- 
ture with  conventional  processing  (uniform  weights,  top  of  Figure  5:  Dolph 
weights,  cop  of  Figure  6).  The  upshot  of  this  is  that  one  buys  improved  adap- 
tive performance  by  filling  in  a sparse  array  at  the  cost  of  additional  elements 

and  processing  load 

A situation  for  which  the  result  is  harder  to  anticipate  is  one  in  which 
a fixed  number  of  elements  are  available  and  the  question  is  whether  to 
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ANGLE  (DEGREES) 

Figure  8.  Comparison  between  the  Adaptive  Performance  of  a 
Sparse  Array  (top)  and  a Filled  Array  of  the  Same 
Ape  rture  (bottom). 


28 


distrihuto  these  to  form  a small  aperture  filled  array  or  a larger  aperture 
sparse  array.  This  case  is  illustrated  in  Figure  9 where  the  adaptive  per- 
formance of  a nine  element  filled  array  (top)  is  contrasted  with  a nine  element 
sparse  array  (bottom).  The  elements  of  the  filled  array  are  spaced  by  a half 
wavelength,  making  the  overall  aperture  of  the  array  four  wavelengths  wide, 
lhe  first  three  elements  of  the  sparse  array  are  spaced  bv  a half  wavelength 
but  the  remaining  six  elements  are  spaced  by  one  and  a half  wavelengths, 
making  the  overall  aperture  of  the  sparse  array  10  wavelengths.  This  sparse 
array  has  all  the  inter -element  spacings  present  in  a filled  10  wavelength  array. 
The  environment  used  for  these  simulations  is  the  same  as  that  in  Figure  3. 

The  result  is  that  tin  effective  angular  resolution  of  the  sparse  array  is  better 
and  the  noise  floor  away  from  discrete  sources  is  somewhat  poorer.  The 
effective  angular  resolution  of  the  sparse  array  is  better  than  that  of  the  filled 
array  by  roughly  the  ratio  of  the  array  apertures  (e.  g.  , the  discrete  sources 
appear  2 to  3 times  as  wide  35  dB  down  for  the  filled  array  than  for  the  sparse 
array).  The  noise  floor  for  the  sparse  array  averages  a few  dB  higher  than 

that  for  the  filled  array,  however. 

The  results  illustrated  to  tlvs  point  have  been  for  small  one -dimensional 
arrays  in  order  to  make  the  displays  simple  and  the  results  easy  to  grasp.  Sim- 
ilar results  have  been  obtained  for  larger  two-dimensionr  1 arrays  of  the  type 

intended  for  data  collection  in  this  program. 

The  cross  is  a two-dimensional  example  of  a low-redundancy  sparsc- 
array:  a filled  array  of  the  same  aperture  would  require  rn  additional  240 
elements,  receivers,  and  processing  channels  and  economically  would  be  out 
of  th  > question.  A 32  element  filled  array  in  two-dimensions  would  have  a 
third  the  angular  resolution  of  the  cross  in  both  azimuth  and  elevation  and 
would  offer  no  important  advantage  in  noise  floor  reduction.  Figure  10  is  a 
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Figure  9.  Comparison  between  the  Adaptive  Perfc  mance  of  a 
Sparse  Array  (bottom)  and  a Filled  Anay  of  the 
Same  Number  of  Elements  (top). 
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Figure  10.  Simulated  Conventional  (top)  and  Adaptive  (bottom)  Sky  Maps 
for  a 32  Flement  Cross  Array.  The  Fnvironment  Consists  of 
Two  Discrete  Source  s of  Unequal  Strength  and  an  Omni-Di- 
rectional Noise  BarKground. 


simulation  for  the  32  element  cross  comparing  the  performance  of  this  array 
with  convi  ntional  beamforming  (top)  and  with  adaptive  beamforming  (bottom). 
Hu  s.  maps  an  formed  by  computing  the  beamfornur  out  nut  for  a 20  by  20 
grid  ol  discrete  azimuth  - elevation  iool<  directions  then  drawing  contours  of 
constant  beamformer  output  through  the  results.  The  environment  consisted  of 
a discrete  "low  angle"  signal  of  0 dB  intensity  (lower  cross  on  figure),  a dis- 
crete "high  angle"  jignal  of  -6  dB  intensity  (upper  cross  on  figure),  and  an 
omni-directional  nois<-  background  producing  a signal  of  -5  dB  intensity  at 
* at  h element.  The  display  is  a radio  map  of  the  sky  in  which  strong  signals 
are  shown  by  double  cross-hatching  (signals  within  3 dB  of  the  strongest 
observed)  and  weaker  signals  by  single  cross-hatching  (signals  3 to  10  dB 
below  t lie  strongest  observed).  With  conventional  processing  the  low  angle 
signal  appears  as  r broad  circular  region  about  4 degrees  in  diameter  (the  3 
dB  boa mwidth  of  the  main  lobe),  surrounded  by  a diffuse  sidetobe  region 
(single  cross-hatching).  The  existence  of  the  high  angle*  signal  is  masked 
because  when  the  beam  is  pointed  in  that  direction,  the  sid°lobes  of  the  patt- 
ern fall  on  the  stronger  source.  In  the  adaptive  map,  the  position  of  the 
strong  low  angle  source  is  much  better  localized  (diametej  of  the  3 dB  con- 
tour is  about  a sixth  as  large),  sidelobe  responses  no  longer  appear,  and  the 
weak  high  angle  source  can  be  clearly  seen. 

These  simulations  suggest  that  sparse  arrays  arc  well  suited  for  use 
with  adaptive  processors,  particularly  in  situations  in  which  the  objective  is 
not  simply  to  detect  a signal  (maximize  SNR)  but  also  to  tell  where  it  is  com- 
ing from  (c.  g.  , CTH  radar  or  passive  HF  direction-finding  applications). 
There  are  an  infinity  of  ways  in  which  a filled  array  can  be  made  sparse:  we 
have  chosen  to  work  with  a particular  class  of  sparse  arrays  (low  redundancy) 
in  which. the  array  is  thinned  to  the  point  where  further  thinning  would 
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, UmiMU-  one  or  more  inter-element  spacing.  present  in  the  parent  lilled 
array.  The  adv.nt.fc  of  this  particular  class  of  sparse  arrays  for  the  pro- 


gram  is  that  the  option  is  prese 


rved  for  producing  high  resolution  maps  ot  the 


environment  by  w 


ell  established  deterministic  beamforming  methods  (Fourier 


transform  of  the  weighted  autocor 


relation  function  for  the  array)  for  compari- 


son with  maps  of  the  same  envir 


onment  produced  by  adaptive  beamforming. 
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4.0  FIELD  INSTALLATION  AND  DATA  COLLFCTION 


A key  objective  of  the  program  has  been  to  collect  data  in  the  real  HF 
environment  and  use  it  to  evaluate  the  benefits  of  adaptive  processing  over 
more  c onventional  techniques.  The  data  for  this  purpose  were  collected  at  an 
established  DARPA  field  site  near  Hudson,  Colorado.  Land,  power,  shelter, 
nrid  much  of  the  electronics  required  were  already  available  at  this  facility  as 
* result  of  other  closely  related  programs.  The  additions  and  modifications 
made  to  this  facility  for  collecting  a data  base  for  adaptive  processing  are  de- 
sc  ribed  below. 

4.1  FACILITY  MODIFICATIONS 

4.1.1  Vertical  Mono-pole  Flements 

A large  circular  array  (3?  elements,  hOO  meter  diameter)  was  available 
at  the  Hudson  site  and  consideration  was  given  to  using  thir  array  without  modi« 
fication.  This  proved  not  to  be  an  attractive  option  because  the  elements  of  the 
circular  array  were  horizontal  crossed  dipoles  designed  for  receiving  signals 
from  the  overhead  ionosphere.  These  elements  were  noc  effective  at  the  low 
elevation  angles  of  practical  interest  in  most  potential  applications  of  adaptive 
processing  (e.g.,  over-the-horizon  radar).  The  principal  requirement  of  the 
data  base  we  were  to  collect  was  that  it  be  representative  of  the  HF  environ- 
ment under  conditions  of  practical  interest  and  consequently  it  was  decided  to 
build  32  new  elements  better  suited  to  the  task.  A vertical  monopole  element 
design  was  available  which  had  been  used  with  success  at  low  elevation  angles 
in  other  applications.  The  radiating  element  is  electrically  short  (5  feet)  to 
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make  it  broadband  and  thick  (1  foot)  to  reduce  the  base  impedance.  A hign 
impedance  base  amplifier  provides  gain  (18  dB)  and  impedance  matching  to 
the  transmission  line  so  as  to  maintain  an  acceptable  system  noise  figure 
over  the  HF  band  <3-30  MHz).  A set  of  32  of  these  elements  were  built  and 
installed  within  the  existing  circular  array.  A sketch  of  the  installation  and 
the  array  elements  is  shown  in  Figure  11.  The  bottom  of  tne  vertical  radiating 
element  was  mounted  about  three  feet  above  ground  level  so  that  the  base 
amplifier  could  be  inspected  with  snow  on  the  ground. 

4.1.2  Cross  Array  Geometry 

The  cross  array  is  a familiar  example  of  a low- redundancy  two-dimen- 
sional array.  It  has  the  virtue  that  the  elements  are  equally  spaced  along 
orthogonal  axes  (in  contrast  to  a circular  array)  which  simplifies  the  process- 
ing software.  It  also  has  the  property  that  the  data  from  the  elements  in  either 
leg  can  be  used  alone  to  study  the  performance  of  a one-dimensional  array  of 
the  same  aperture  as  the  parent  two-dimensional  array.  For  these  reasons 
the  cross  array  was  chosen  for  use  in  the  program.  A sketch  of  the  installed 
array  is  shown  in  Figure  11  and  a dimensional  plan  view  is  shown  in  Figure  12. 
The  inter-element  spacing  is  20  meters  (a  half-wavelength  at  7.  S MHz).  The 
east-west  leg  of  the  array  is  16  elements  long  (a  300  meter  aperture)  and  the 
north-south  leg  is  17  elements  long  ( a 320  meter  aperture'. 

4.1.3  Receiving  and  Data  Acquisition  System 

The  receiving  and  data  acquisition  system  used  to  'ollect  the  adaptive 
data  base  is  shown  in  Figure  13.  The  new  instrumentation  required  to  assem- 
ble this  system  from  the  equipment  already  available  at  the  Hudson  site  was 
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*—  RING  RECEIVING  ARRAY 
(32  HORIZONTAL  CROSSED  DIPOLES) 


Figure  12.  A Plan  View  oC  the  Hudson  Field  Site  Showing  the 

Orientation  and  Inter-Element  Spacing  of  the  Cross 
A r ray. 
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a sot  of  32  "Q"  receiver  modules,  the  associated  power  dividers,  and  the  HF  ) 
distribution  circuicry.  These  receivers  are  identical  to  the  original  set  of  32 
but  the  HFO  signal  is  shifted  by  90*.  A pair  of  receivers  is  used  for  each 
element,  providing  the  complex  samples  required  for  adaptive  processing. 

The  64  receiver  channels  are  all  sampled  simultaneously  and  the  results  held 
while  an  A/I)  converter  steps  from  channel  to  channel  and  transfers  the  data 
to  the  minicomputer  memory  (dish  storage).  The  process  is  repeated  peri- 
odically until  the  computer  memory  is  full:  data  collection  then  stops  while 
the  stored  information  is  recorded  on  digital  tape  (four  minutes  required  to 
accomplish  this).  The  memory  will  hold  90112  words  which  corresponds  to  22 
blocks  of  data,  each  block  containing  64  time  samples  of  the  entire  array,  each 
sample  being  a set  of  64  receiver  channel  voltages.  The  data  base  for  this  pro- 
gram was  collected  r.t  the  rate  of  hO  samples  of  the  cnti-e  array  per  second. 
This  corresponds  to  one  data  block  every  1.28  seconds  and  the  22  block  capa- 
city of  the  memory  in  28.16  seconds.  The  data  acquisition  sequence  thus  al- 
ternated between  half-minute  periods  of  data  collection  and  four  minute  periods 
of  data  recording. 

4.2  TEST  GFOMETR  Y 

Pata  were  collected  with  and  without  cooperative  C W transmissions  from 
an  SRI  field  site  at  Tost  Hills,  California.  The  SRI  transmissions  provided 
signals  under  realistic  but  controlled  conditions  and  it  is  only  data  from  this 
portion  of  the  data  base  which  have  been  analyzed  in  this  report.  The  geometry 
of  the  path  from  Lost  Hills,  California  to  Hudson,  Colorado  is  shown  in  F igure 
14.  The  length  of  the  path  is  1433  km  and  the  azimuth  (measured  from  Hudson) 
is  2 ^4 . 5 degrees . 
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Figure  14.  Geometry  of  Cooperative  Tests  with  SKI  Providing 
a Controlled  Interfering  Signal. 
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4.3  DATA  COLLECTION 

J 

The  purpose  of  the  data  collection  effort  was  to  obtain  a samplinR  of  the 
HT  environment  over  a ran^e  of  frequencies  and  interference  conditions. 

• / 

Tahli  1 summarizes  the  environmental  data  on  tape. 


TABLE  I:  ENVIRONMENTAL  DATA  BASF 


Pate 

I req. 

F nvironment 

Records 

18  Nov.  1973 

7 MHz 

Amateur  Radio  Band  (heavy, 
multiple  source  interference) 

44 

19  Nov.  1973 

10  MHz 

Cooperative  SRI  Transmissions 
(multiple  hop,  element  SNR  ->1) 

17  f. 

15  MHz 

Cooperative  SRI  Transmissions 
(single  hop,  element  SNR  > i) 

176 

IS  MHz 

Cooperative  SRI  Transmissions 
(single  hop,  element  SNR  < 1) 

176 

28  .Inn.  1974 

15  MHz 

International  Broadcast  Band 
(heavy,  multiple  hop  interfer- 
ence) 

924 

Each  data  record  contains  64  complex  samples  of  each  array 

element  and 

represents 

1.28  seconds 

in  time. 

43 


/ 


S.o  ANALYSIS  OF  EXPERIMENTAL  DATA 


r1#  1 s If.  N A I , FN  VI HONMFNT 

The  data  analyzed  in  this  section  was  taken  as  representative  of  the 
least  complex  adaptive  environment  likely  to  be  encountered  in  practice,  a 
strong  C'VV  interferer  arriving  via  a single  hop  path.  The  interfering  C VV  sig- 
nal was  provided  by  SRI  (see  Figure  14  for  test  geometry)  at  a frequency  o* 

1 .720  MHz.  This  frequency  was  chosen  so  as  to  lie  just  below  the  maximum 
usable  frequency  (MUF)  for  the  path  at  the  start  of  the  data  taken  to  insure 
that  only  one-hop  signals  would  be  observed.  The  one-hop  signal  will  of 
course  he  composed  ofhigh  and  low  angle,  ordinary  and  extraordinary  modes 
hut  by  operating  near  the  MUF  we  hoped  to  minimize  their  angle  of  arrival 
differences  and  obtain  an  incident  signal  on  the  array  who  h approximated  a 
single  point  source  in  the  far  field.  Such  a source  would  produce  a plane  wave 
incident  on  the  array,  a common  assumption  in  numerical  simulation  studies 
of  adaptive  processors  and  therefore  an  important  situation  to  study  experi- 
mentally. Ionospheric  conditions  were  such  that  the  SRI  signal  arrived  at  the 
array  at  an  elevation  angle  of  about  IS  degrees  and  at  an  azimuth  angle  of  about 
2S4  degrees.  The  SNR  in  the  full  17  kHz  bandwidth  of  the  receiver  was  of 
order  10  dB  (measured  at  the  element,  prior  to  array  gain). 


S.2  ONE-DIMENSIONAL  RESULTS 

The  simulations  presented  in  Section  3.0  emphasized  one -dimensional 
arrays  because  the  results  from  such  arrys  are  more  convenient  to  display 
in  detail  and  because  the  concepts  involved  are  qualitatively  the  same  in 
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pithrr  one  dimension  or  two.  One  of  the  advantages  of  collecting  data  on  a 
cross  array  is  that  the  outputs  from  one  leg  at  a time  can  he  processed  after 
the  fact  to  give  results  for  a one-dimensional  array  and  this  has  been  done  for 
comparison  with  one -dimensional  simulations.  Since  the  SHI  signals  arrive 
from  the  west,  the  north-south  leg  of  the  array  was  used  to  form  a 17  element 
broadside  array  and  subsets  of  this  array  were  proc  essed  for  comparison  with 
the  simulations.  The  results  for  a 1.28  second  data  block  are  shewn  in  F igure 
l . The  plot  at  the  top  of  this  figure  shows  the  power  output  of  the  adaptive 
signal.  This  result  is  for  an  array  of  IS  elements  spaced  20  meters  apart. 

The  3 dB  conventional  beamwidth  of  the  array  (at  IS.  729  MHz)  is  4 degrees 
and  this  angular  width  is  indicated  on  the  figure  at  the  3 db  down  level.  As 
anticipated  from  the  simulations  (Figures  S and  6 of  Sec  tion  3.0),  the  defini- 
tion of  the  position  of  the  discrete  source  is  much  sharper  with  adaptive  pro- 
cessing than  with  conventional  (by  a factor  of  about  1 8 at  the  3 HB  down  level 
in  this  example).  The  side  lobe  response's  are  well  below  the  13  dB  down  level 
c'xpc'ctc'd  for  conventional  processing.  The  plot  at  the  bottom  of  this  figure 
shows  the  result  of  removing  9 of  the  original  18  elements  of  the  one-dimen- 
sional array  to  form  a sparse  array  of  nearly  the  same  aperture  but  less  than 
half  the  number  of  elements.  As  anticipated  from  the  simulations  (Figure  7 
of  Section  3.0),  the  angular  width  of  the  discrete  source  is  essentially  the 
same  for  the  sparse  array  as  for  the  filled  array.  The  conclusion  is  that  the 
simulations  have  proven  to  be  in  satisfactory  qualitative  agreement  with  re- 
sults obtained  in  the  teal  HF  environment. 
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AZIMUTH  ANSLE(DEGREES) 


Figure  IS.  Comparison  of  Azimuth  Maps  of  SRI  Signal  (IS.  729  MHz) 
with  Filled  Array  (top)  and  Sparse  Array  (bottom).  The 
Sparse  Array  has  fewer  elements  hut  Comparable  Aper- 


5.3  TWO-DIMENSIONAL  RESULTS 


An  attractive  feature  of  a two-dimensional  array  for  a study  of  adaptive 
processing  performance  is  that  it  permits  one  to  map  the  radio  sky  to  show 
where  the  interfering  signals  are  coming  from  and  then  to  show  how  the  noise 
floor  reduction  varies  as  a function  of  position  relative  to  these  signals.  The 
noise  floor  reduction  is  thus  not  presented  as  a single  number  but  rather  as  a 
contour  map  showing  the  improvement  in  the  noise  floor  as  a function  of 
azimuth  and  elevation  in  the  sky.  We  have  called  this  display  an  'improvement" 
map.  The  improvement  in  the  noise  floor  achieved  in  a given  direction  is  the 
difference  between  the  power  output  of  the  conventional  processor  and  power 
output  of  the  adaptive  processor.  The  improvement  map  is  a contour  plot 
showing  how  the  noise  floor  reduction  varies  as  a function  of  azimuth  and 
elevation  across  the  radio  sky.  Since  the  conventional,  the  adaptive,  and  the 
improvement  map  are  closely  related  we  have  displayed  them  together  as  a 
set  for  each  data  block,  along  with  a plot  of  the  frequency  spectrum  of  the  re- 
ceived signal  from  which  the  element  SNR  for  that  data  block  can  be  read. 

The  SRI  signal  was  characterized  by  a fading  period  of  about  5 seconds: 
a 1.28  second  block  length  was  short  enough  so  that  the  signal  level  did  not 
change  significantly  over  that  interval.  In  order  to  examine  conditions  over 
a range  of  signal  levels  a series  of  11  contiguous  data  blocks  were  processed, 
corresponding  to  about  14  seconds  in  time  and  covering  about  3 fading  periods. 

Results  for  the  first  data  block  are  shnwti  in  Figure  16.  The  azimuth 
and  elevation  scales  for  the  three  maps  are  the  same  and  show  the  sector  of 
the  sky  containing  the  interfer.  The  scales  are  linear  in  "sine-space",  mak- 
ing them  non-linear  in  angle  in  such  a way  that  the  foreshortening  of  the 
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Ficurc  16.  Conventional  Map.  Adaptive  Map,  Improvement  Map. 

and  Frequency  Spectrum  for  Data  Block  Number  1. 
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aperture  which  occurs  as  the  array  is  steered  off-boresight  (zenith)  is  com- 
pensated for  on  the  display.  The  effect  is  to  keep  the  shape  of  the  antenna 
pattern  constant  on  the  display  regardless  of  position  and  consequently  make 
point  sources  at  low  elevation  angles  appear  as  round  images  rather  than  as 
elliptical  images  extended  along  the  elevation  angle  axis.  The  conventional 
map  shows  the  SRI  signal  at  about  1 S degrees  elevation  angle  and  about  2 r>4 
degrees  azimuth.  Gray  scales  are  used  to  show  signal  level.  The  darkest 
irca  represents  the  first  three  dB  down  from  the  peak  signal  and  shows  that 
tin  t onventional  half-power  beamwidth  of  the  arra^  is  about  5 degrees.  The 
theoretical  value  for  a true  point  source  is  4 degrees:  the  slightly  greater 
measured  beamwidth  is  possibly  due  to  the  finite  angular  dimensions  of  the 
source.  The  peak  sidelobe  responses  of  the  array  are  strongest  along  two 
orthogonal  arms  which  mirror  the  geometry  of  the  cross:  the  peak  sidelobe 
level  is  of  order  5 dB  below  the  main  beam  response.  Thus  relatively  high 
sidelobe  level  is  characteristic  of  sparse  arrays  with  conventional  processing 
and  would  not  lie  helped  by  applying  a fixed  set  of  weights  (e.  g.  , Dolph- 
Tchebyscheff)  in  the  beamforming  process.  The  sidelober.  could  be  reduced 
by  filling  the  array  'requiring  an  additional  240  elements  and  receivers)  and 
then  applying  a fixed  set  of  weights  but  this  would  not  improve  the  angular  re- 
solution. The  adaptive  map  does  not  exhibit  the  high  sidelobe  responses  of 
the  conventional  map  and  the  "main  beam"  signal  definition  is  somewhat  better. 
It  is  probable  that  the  signal  at  this  time  actually  consists  of  two  or  more  modes 
of  comparable  strength  separated  in  angle  (or  perhaps  a single  mode  'Scintill- 
ating" in  angle  during  the  1.28  seconds  used  to  collect  the  data).  The  contours 
for  the  adaptive  map  have  been  drawn  in  terms  of  dB  below  the  power  level 
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obtained  when  the  adaptive  processor  is  steered  in  the  direction  of  the  inter- 
feres This  makes  possible  a direct  comparison  between  the  shapes  of  con- 
ventional and  adaptive  sky  map  patterns  but  it  does  obscure  the  fact  tha>  the 
maximum  power  output  of  the  adaptive  processor  was  18  dB  less  than  the 
maximum  power  output  of  the  conventional  processor  in  the  source  direction. 
f.  or  an  ideal  point  source  the  two  processor  outputs  should  be  the  same  and 
tentatively  attribute  the  discrepancy  to  the  finite  dimensions  of  the  source: 
an  equivalent  statement  is  that  the  signal  incident  on  the  array  departs  siRuifi- 
c-rtntly  from  a plane  wave.  The  result  is  a lowering  of  the  noise  (interference) 
flooi  with  the  adaptive  processor  when  steered  in  the  general  direction  of  the 
interfere,  in  contrast  to  simulation  results  based  on  ideal  plan.-  wave  signals 
which  predict  no  noise  floor  reduction  when  looking  at  the  interfered  The 
improvement  map  itself  shows  noise  floor  reductions  between  20  and  30  dB 
over  much  of  the  sky,  though  there  are  regions  where  little  improvement  is 
obtained.  The  SNR  is  an  important  parameter  in  the  comparison  of  field  data 
with  simulations:  we  have  estimated  this  quantity  from  the  power  spectrum 
of  the  signal  using  one  receiver  pair  (I  and  Q channels  for  one  array  element). 
The  signal  energy  is  confined  to  within  a few  Hertz  of  the  center  frequency. 
The  frequency  bins  outside  this  region  have  been  taken  as  noise  and  used  to 
compute  an  average  noise  level.  It  has  then  been  assumed  that  this  noise 
lev-1  is  also  representative  of  the  noise  level  within  the  frequency  bins  con- 
taining signal  (horizontal  dashed  line  on  the  figure  separating  signal  from 
noise).  The  total  noise  power  and  the  total  signal  power  across  the  spectrum 
have  then  been  computed  and  these  numbers  used  to  compute  the  integrated 
SNR  shown  on  the  figure  (12  dB).  AH  of  the  noise  power  in  the  17  kHz  band- 
pass the  receiver  is  represented  in  the  displayed  spectrum  because  of 
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aliasing  into  the  ±25  Hz  window  displayed  on  this  figure.  The  12  dB  SNR  for 
data  block  § 1 is  one  of  the  highest  of  this  series,  even  though  the  adaptive  and 
the  improvement  maps  are  among  the  least  impressive.  This  leads  us  to  be- 
lieve that  variations  in  mode  structure  from  map  to  map  are  playing  an  impor- 
tant role  in  the  performance  of  the  adaptive  processor. 

Data  block  /2  is  shown  in  Figure  17  and  is  by  all  measures  the  least  im- 
pressive examples  of  adaptive  (or  conventional)  processing.  The  element  SN'R 
is  the  lowest  observed  (8  dR),  the  conventional  map  is  poorly  defined  and  has 
an  uneven  sidelobe  structure,  the  adaptive  map  shows  somewhat  better  signal 
structure  definition  and  freedom  from  sidelobe  responses,  and  the  imorove- 
ment  map  shows  only  modest  noise  floor  reductions  (10  - 20  dB  across  most 
of  the  sky). 

Data  block  / 3 (Figure  18)  exhibits  the  highest  SNR  observed  in  this  series 
(13  dB),  more  regular  sidelobe  structure  in  the  conventional  map,  the  pre- 
sence of  a more  compact  signal  structure  on  the  adaptive  map,  and  noise  floor 
reductions  in  the  20  to  30  dB  range  on  the  improvement  map. 

Data  block  / 4 (Figure  19)  is  one  of  the  best  examples  obtained  of  the  bene- 
fits of  adaptive  processing  at  HF.  The  conventional  map  is  regular,  the  adap- 
tive map  provides  a much  sharper  picture  of  the  interfering  signal,  and  the 
improvement  map  shows  noise  floor  reductions  in  the  30  to  40  ciB  range.  The 
element  SNR  for  this  data  block  was  not  outstanding  (10  dB). 

Data  block  / f Figure  20)  is  the  best  example  of  the  sort  of  results  which 
would  have  been  expected  from  simulations  with  a single  plane  wave  inter- 
ferer.  The  conventional  map  exhibits  the  classical  cross  sidelobe  pattern, 
the  adaptive  map  shows  a single,  sharply  defined  source,  and  the  improve- 
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Figure  18.  Conventional  Map,  Adaptive  Map,  Improvement  Map, 
and  Frequency  Spectrum  for  Data  Block  Number  3. 
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19.  Conventional  Map,  Adaptive  Map, 
• and  Frequency  Spectrum  for  Data 

Improvement  Map, 
Block  Number  4. 
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Figure  20.  Conventional  Map,  Adaptive  Map,  Improvement  Map, 
and  Frequency  Spectrum  for  Data  Block  Number  5. 
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menl  map  shows  noise  floor  reductions  of  order  30  - 40  dll  over  much  of  the 
sky.  The  3 dll  width  of  the  "main  lobe"  response  on  the  adaptive  map  ts 
about  6 times  sharper  than  on  the  conventional  map.  This  indicates  that  the 
interferer  approximated  a point  source  to  within  a degree  of  angular  spread: 
the  improvement  map  also  shows  the  dip  in  improvement  expected  from  sirm- 
iations  when  the  adaptive  processor  is  steered  in  the  direction  of  a point  source 
interferer.  Tver,  here,  however,  the  output  of  theadaptive  processor  is  more 
than  20  dll  below  that  of  the  conventional  processor  when  pointed  in  the  nominal 
direction  of  the  ir.terlerer  whereas  simulations  with  an  ideal  point  source  >n- 
terferer  would  predict  equal  output  in  both  cases.  The  element  SNR  for  this 
data  block  is  11  dli  and  thus  is  intermediate  between  the  highest  and  lowest 

observed  in  this  series. 

Data  block  f 6 through  / 11  (Figures  21  through  26)  cornolete  the  series  and 
are  presented  to  show  that  while  the  details  vary  from  block  to  block,  the  first 
S data  blocks  are  representative  of  what  can  be  expected  under  these  conditions. 
These  later  data  block,  are  self-explanatory  and  no  further  discussion  will  be 

given. 

To  summarize,  the  preceding  data  set  has  she  vn  that  adaptive  processing 
of  two-dimensional  HF  arrays  consistently  produces  improved  definition  of  the 
angle -of -arrival,  freedom  from  sidelobe  responses,  and  significant  noise  floor 
reduction  compared  to  conventional  processing.  The  absolute  magnitudes  of 
the  improvement,  ere  of  course  dependent  upon  the  particular  antenna  array 
used  and  the  particular  set  of  field  data  chosen  for  analysis.  The  importance 
of  what  has  been  ,-ccomplished  is  that  it  has  been  shown  that  simulations  using 
existing  algorithms  are  satisfactory  in  most  respects  fot  predicting  the  per- 
formance of  two-dimensional  adaptive  processors  in  the  real  HF  environment. 
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Figure  21  . Conventional  Map,  Adaptive  Map,  Improvement  Map 
ard  Frequency  Spectrum  for  Data  Block  Number  6. 
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Figure  22.  Conventional  Map,  Adaptive  Map,  Improvement  Map, 
and  Frequency  Spectrum  for  Data  Block  Number  7. 
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F ii»ure  2 3.  C onvontional  Map,  Adaptive  Map,  Improvement  Map, 
and  Frequency  Spectrum  for  Data  Flock  Number  8. 
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Figure  2 6.  C conventional  Map,  Adaptive  Map,  Imp-ovement  Map, 
and  Frequency  Spectrum  for  Data  Block  Number  10. 
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Figure  2(>.  Conventional  Map,  Adaptive  Map,  Improvement  Map 
and  Frequency  Spectrum  for  T)ata  Block  Number  !!• 
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The  one  aspect  of  the  results  which  differ  markedly  from  that  expected  is  the 
fact  that  the  adaptive  processor  output  is  consistenly  we'l  below  (of  order  20 
(IB)  Ihe  output  of  the  conventional  processor  when  both  are  steered  in  the  di- 
rection of  the  interferer,  an  effect  tentatively  attributed  (o  Ihe  finite  angular 
extent  of  the  incident  signal,  and  the  greater  sensitivity  of  adaptive  processors 
to  such  departures  ‘liar,  of  conventional  processors. 

The  last  topic  to  be  discussed  in  this  section  is  the  stability  of  the  medium, 
as  it  affects  the  noise  floor  reduction  achieved  through  adaptive  processing. 
The  cost  ot  implementing  an  adaptive  processor  in  real  time  depends,  among 
other  things,  upon  how  frequently  one  must  update  the  adaptive  weights.  The 
impru  oment  maps  shown  to  this  point  were  computed  with  weights  optimized 
for  each  data  block.  If  instead,  the  data  for  each  block  were  analyzed  with 
•he  weights  computed  for  the  preceding  block,  the  noise  floor  reduction  would 
ue  less:  the  question  is  "how  much  less?".  Figures  27  through  36  compare 
the  improvement  maps  obtained  using  current  weights  (left)  and  using  weights 
1 block  old  (right!.  The  results  vary  somewhat  from  block  to  block  and  from 
one  position  in  the  map  to  another  but  in  round  numbers  there  is  a 10  dB  loss 
in  improvement  if  weights  for  the  preceeding  1 second  data  block  are  used. 

By  way  of  comparison,  Griffiths  (reference  6)  found  using  8 second  data  blocks 
that  the  penalty  paid  for  using  weights  optlmlzedfor  the  preceding  block  was 
also  of  order  of  10  dB.  It  must  be  cautioned  that  the  experimental  conditions 
were  quite  different  in  these  two  sets  of  measurements  and  that  the  data  base 
in  both  cases  was  very  limited.  However,  it  appears  that  shortening  the  block 
length  from  8 to  1.28  seconds  did  not  reduce  the  SNR  penalty  paid  for  using 
weights  optimized  for  the  preceding  block  (10  dB).  It  is  clearly  important 
to  update  the  weights  for  each  block  but  it  remains  to  be  seen  how  the  average 
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igure  36.  Improvement  Maps  for  Data  Block  II  with  Optimum  Weights  deft)  and 

rig! 


SNR  for  a block  depends  upon  block  length.  The  processing  load  will  decrease 
as  the  block  length  increases  (the  time  consuming  operation  is  the  matrix  in- 
version required  once  per  block  to  update  the  weights).  The  data  on  tape  should 
be  reprocessed  with  progressively  longer  block  lengths  to  see  how  the  SNR 
improvement  varies  with  block  length. 

There  are  several  other  areas  of  the  analysis  which  remain  unexplored. 

The  data  presented  in  this  section  should  be  reprocessed  after  first  filtering 
in  the  frequency  domain  to  improve  the  SNR  (the  improvement  expected  is  about 
lo  dB).  OTH  radar  systems  will  have  to  do  this  to  achieve  target  detectability 
and  it  is  of  interest  to  see  what  happens  to  the  improvement  stability  times  as 
the  bandwidth  becomes  less  than  the  signal  bandwidth.  Increased  SNR  will  al- 
so permit  us  to  make  a more  stringent  test  of  the  nulling  capability  of  a prac- 
tical two-dimensional  array,  in  that  sidelobe  responses  are  now  undetectable 
in  the  adapted  maps. 

The  data  taken  with  both  1 and  2 hop  signals  from  SRI  presented  simultane- 
ously should  be  analyzed  to  examine  the  performance  of  an  adaptive  processor 
under  the  common  situation  in  which  a strong  interferer  arrives  via  two  or 
more  modes  well  separated  in  elevation  angle.  The  two-hop  signal  could  be 
viewed  as  a desired  signal  and  the  one-hop  as  the  interferer,  for  example, 
a two-dimensional  array  should  be  able  to  deal  with  this  situation  effectively 
whereas  a one-dimensional  array  (azimuth  only)  could  no!. 

Finally,  the  data  taken  under  conditions  of  heavy,  multiple  source  inter- 
ference (Ham  bands,  international  broadcast  bands)  should  be  used  to  examine 
the  processor  performance  under  conditions  where  the  number  of  sources  are 
greater  than  the  number  of  array  elements  (which  determines  the  number  of 
degrees  of  freedom  available  and  the  number  of  nulls  which  can  be  steered). 
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6.0  SUMMARY  AND  CONCLUSIONS 


Hu  DARPA  phased  array  facility  near  Hudson,  Colorado  lias  ueen 
modified  to  permit  the  collection  of  a data  base  on  a two-dimensional  HP  array 
in  a format  suitable  for  adaptive  processing.  This  has  included  the  fabrication 
and  installation  of  32  receivers,  32  vertical  monopole  elements  with  base  ampli 
fiers,  and  the  development  of  software  for  controlling  the  data  acquisition  and 
recording  system.  Data  have  been  collected  under  a variety  of  conditions  and 
in  a form  such  that  adaptive  processing  can  be  done  in  azimuth,  in  elevation, 
in  frequency,  or  in  any  combination  of  these.  Analysis  software  has  been 
written  to  reduce  the  data  with  conventional  beamforming  algorithms  and  with 
existing  adaptive  beamforming  algorithms  (direct  matrix  inversion).  Analysis 
software  has  also  been  written  to  produce  simulate  1 data  tapes  for  model 
environments  and  model  array  configurations:  these  have  been  used  to  verify 

the  data  reduction  software  and  to  study  the  < ffect  of  array  configuration  on 
adaptive  antenna  performance.  The  field  data  tapes  have  been  made  available 
to  the  DARPA  community  and  a limited  amount  of  analysis  of  these  tapes  has 
been  done  at  Raytheon.  The  principal  findings  of  this  analysis  are  the  follow- 
ing: 

1.  Verification  that  the  11F  environment  will  in  fact  permit  several 
of  the  b<-n"fits  anticipated  from  numerical  simulations;  notably, 
a.  Noise  Floor  Reduction:  reduction  of  the  noise  floor  by 

30  dB  or  more  over  much  of  the  sky.  This  was  achieved 
with  pattern  control  alone  (adapting  in  azimuth  and  eleva- 
tion but  not  in  frequency)  against  an  int<  rferer  of  modest 
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strength  (10  e|  B SNR  nicasuri  'I  at  tin  clement), 
h.  Snli  lulu  Reduction  in  alternate-  wa*  of  viewing  the  noise- 

floor  re-due t ion  achieve-d  is  that  the-  ette-itiv  - sielelobe  levt  Is 
cf  the1  antonna  we  re  reduced  by  more  than  30  dB.  F or  the- 
elata  examined  to  date-  (a  single  dominant  -nte  rferor  with 
SNR  of  10  elB  at  the  element)  the  effect  of  antenna  side  lobe-s 
we  re  not  me-re-ly  reduced  but  -ere  eliminated,  in  that  siele-- 
lobe  re'sponscs  were  undetectable  in  the  adaptively  processed 
data . 

c.  EO-ctive  Beamwidtlr  the-  effective  3 dB  beamwidth  of  the- 
ante  nr.a  array  (after  processing)  was  reduced  by  as  much 
as  a factor  of  6 when  used  to  map  isolated  signals  having 
ele-me  nl  SNR's  of  order  10  dB.  That  is,  the  apparent  anpu- 
lar  width  of  discrete  source  s on  sky  maps  was  much  narrower 
with  i elaptive  processinp  than  with  conv  c-ntiemal  proce-ssing. 

2.  Measurement  of  the  following  properties  of  the  environment  whic  h 
impact  directly  on  the  usefulness  of  adaptive  pie>cessinp  at  IIF . 
a.  Stability  of  the-  Environment:  the  environment  chanpes  rapidly 
eno  ill  so  that  if  the  optimum  adaptive-  weights  are  replaced 
by  weights  1 second  old,  the  loss  in  noise  floor  reduction  is 
order  10  dB.  For  a situation  in  which  the  noise.-  floor  reduc- 
tion with  optimum  weights  was  3l  dB,  this  me-ans  that  the- 
reduction  will  be  only  of  order  20  ciB  if  weights  1 second  aid 
are  used.  The  upshot  of  this  is  to  .uggest  that  in  practice 
the-  adaptive  weights  will  have  to  be  recomputed  for  each  data 
block:  it  remains  to  be  seen  how  long  individual  data  blocks 
can  be  without  significant  SNR  penalty. 
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b.  Limitations  of  Plane-Wave  Idealization:  t*>r  an  ideal  antenna 
(identical  element  patterns,  no  calibration  errors),  ideal 
terrain  (no  scattering  from  local  structures),  and  ideal  iono- 
spheric reflection  (single  discrete  mode),  the  amplitude  of 
the  signal  would  he  uniform  across  the  array  and  the  phase 
would  vary  linearly.  For  this  case  the  output  of  conventional 
and  ac\  ptivo  processors  would  be  the  same  when  steered  to 
the  source  direction.  In  practice,  however,  it  was  found  that 
the  tit  parturt  s from  ideal  were  enough  {e.  g.  , RMS  amplitude 
variations  of  a few  t<  ns  of  percent  across  the  array)  to  cause 
ti  e output  of  the  adaptive  processor  to  he  M to  10  dB  less  than 
the  output  of  the  conventional  processor  when  steered  in  the 
source  direction.  This  result  was  unexpected,  though  in 
retrospect  it  could  have  been  predicted  from  the  RMS  varia- 
tion from  plane  wave  conditions  reported  for  other  large  HF 
a-rays  and  numerical  simulations  of  th<  impact  of  such  varia- 
tions on  adaptive  processor  performance.  It  is  not  known  to 
what  extent  the  observed  departure  from  the  plane  wave 
idealization  is  due  to  residual  system  calibration  errors, 
scattered  fields  from  local  irregularities,  or  irregularities 
introduced  by  the  ionospheric  reflection  process.  The 
recommendation  is  that  future  data  collection  efforts  he  pre- 
ceded by  a test  in  which  a beacon  transmitter  is  flown  across 
the  array  to  verify  the  calibration  of  the  system  as  a whole 
(antenna  and  local  terrain).  It  should  he  emphasized  that  the 
impact  of  departures  from  ideal  plane  wave  conditions  depends 
upon  the  SNR  of  the  discrete  signal  and  may  not  be  a problem 
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tor  1 hr  low  SNR  signals  of  common  interes-  in  applications. 
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3.  The  effort  of  array  configuration  on  adaptive  processor  perform- 
ance was  considered  from  the  viewpoint  of  the  systt  m designer 
seeking  to  place  the  available  elements  to  best  advantage:  low 

redundancy  sparse  arrays  (be  they  one  or  two  c'imensional)  art* 
suggested  as  promising  configurations  for  making  the  most  of  the 
* lemcnts  and  processing  capability  available.  U has  not  been 
shown  rigorously  that  such  arrays  are  "best"  in  any  specific 
sense  but  numerical  simulations  and  comparisons  with  field  data 
mak(  it  plausible  for  situations  in  which  'he  goal  is  not  merely 
to  detect  a signal  (maximize  SNR)  but  also  to  tell  where  it's 
coming  from  (maximize  angular  resolution). 


' 
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